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We present continuum observations of the Perseus Bl-E region from the Herschel Gould Belt Survey. These Herschel data reveal 
a loose grouping of substructures at 160 - 500 fim not seen in previous submillimetre observations. We measure temperature and 
column density from these data and select the nine densest and coolest substructures for follow-up spectral line observations with the 
Green Bank Telescope. We find that the Bl-E clump has a mass of ~ 100 M Q and appears to be gravitationally bound. Furthermore, 
of the nine substructures examined here, one substructure (B1-E2) appears to be itself bound. The substructures are typically less than 
■ a Jeans length from their nearest neighbour and thus, may interact on a timescale of ~ 1 Myr. We propose that Bl-E may be forming 

a first generation of dense cores, which could provide important constraints on the initial conditions of prestellar core formation. Our 
results suggest that Bl-E may be influenced by a strong, localized magnetic field, but further observations are still required. 
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1 . Introduction (i.e., clustered or isolated star formation and low-mass or high- 
mass star formation). Furthermore, theoretical studies also sug- 
Molecular clouds are highly structured regions of dust and gas. gest that additional processes, such as radiation feedback (i.e., 
They contain dense, small-scale, star-forming "cores" (< 0. 1 pc) iKrumholz et al.l2010h or large-scale sho cks associated with con- 
that are usually cl ustered into larger-sca l e clumps and organized verging flows (i.e.. Heitsch et al.ll201 ll) . can influence molecu- 
k>( ! along filaments (IWilliams et"aLTl2000t IWard-Thompson et all lar cloud structure and thus, core formation. Additionally, re- 
■ 120071: IdT Francesco et al.ll2007l: lAndre et al.ll2010t> . The cause of cent Herschel studies have shown that filaments are very promi- 
& \ this hierarchical structure, where the large-scale clouds (~ 10 nent in star forming regions and likely have an important role in 
pc) with moderate densities (~ 10 2 cm -3 ) produce filaments, the formation of dense substructures in mo lecular clouds (e.g., 
clumps, and den se cores at higher densit ies (> 10 4 cm 3 ), is not lAndre et aLll2010t lArzotimanian et al 

well understood (Bergin & Tafalla 2007). c c ... . . . , . 

* b Surveys of core populations in molecular clouds indicate 

Molecular clouds form dense cores when diffuse gas that cores are confin ed t o regions of high column density. 

is compressed. Current theories for core formation pri- iJohnstone et al.1 (120041) and iKirk et all (120061) found a relation- 

marily focus on two ideas: that cores for m via ( 1 ) tur- ship between core occurrence and extinction for the Ophiuchus 

bulent compression of d iffuse gas (e.g ., iLarsonl 1 198 it and Perseus clouds, respectively, suggestin g that there is a core 

iMac Low & Klessenl l2004t iDib et all I2009I) . or (2) the mo- formation threshol d at A v > 5. Similarly. ILada et all (12009-1) and 

tion of neutral material between magnet i c field lines or am- lAndre et"ail J2010) each compared two clouds with different de- 

bipolar diffusion (e.g.. [Mestel & Spitzerf 119561: iMouschoviasI grees of star formation activity and each found that the more 

1976; Kunz & Mouschovias 2009). Both mechanisms and grav- active cloud was composed of higher column density material 

ity likely play important roles in star formation, but it is un- (by a factor of ~ 10) than the quiescent cloud. These studies 

clear whether one mechanism would dominate in all situations emphasize that cores require a minimum column density (ex- 
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tinction) of material to condense from the bulk cloud (see also 
iHeiderman et al.ll2010l) . 

The precursors to cores are difficult to identify. Dense cores 
are often influenced by processes such as nearby outflows 
and radiation feedback from a previous epoch of nearby star 
formation, and observations of them cannot be used to con- 
strain Jhe_d^iarriicpro£erties of the initial core-forming mate- 
rial dCurtis & Richerll201 ll) . Without knowing the initial condi- 
tions and processes that cause cores to form from diffuse gas, 
we cannot accurately model their formation or evolution. Thus, 
identifying and analyzing a core forming region without earlier 
episodes of star formation would be exceedingly useful to con- 
strain how molecular clouds form star-forming substructures. 

In this paper, we use observations from the Herschel Gould 
Belt Survey to explore a ~ 0.1 deg 2 clump r oughly 0.7° east of 
the Bl clump (Bac hiller & Cernicharof l986) within the Perseus 
molecular cloud. This region, hereafter called B 1 -E, has high ex- 
tinction {A v > 5) similar to the core formation threshold. Despite 
this high extinction, previous submillimetre and infrared contin- 
uum observations suggest that Bl -E contains neither dense cores 
nor young stellar objects (e.g.. Enoch et all 120061: iKirk et all 
2006: Ifergensen et alJl2007t lEvans et al.ll2009l) . In contrast, our 
Herschel observations show substructure that was not detected 
by these other far-infrared and submillimetre continuum studies. 
Furthermore, we use recent Green Bank Telescope (GBT) obser- 
vations to quantify in part the kinematic motions of the densest 
substructures seen in the Herschel data. 

We propose that Bl-E is forming a first generation of dense 
cores in a pristine environment. In Section [2] we describe our 
Herschel and GBT observations. In Section [3] we describe the 
properties we derive from our data. In Section [4] we discuss the 
implication of our results and compare our observations to theo- 
retical models. Finally, in Section[5]we summarize the paper. 

2. Data 

2.1. Herschel Observations 

The western half of Perseus, including Bl-E, was observed 
by the Photodetector Array Camera and Spectrometer (PACS; 
Pogli tsch et al.ll20 10|) and the Spectral and Photometric Imaging 
Receiver (SPIRE: iGriffin et al.ll2010h as part of the Hersche l 
Gould Belt Survey (An dre & Saracenol2005tlAndre et al.l2010h . 
The Herschel observations were taken in "parallel mode", result- 
ing in simultaneous coverage at 5 wavelength bands, with the 70 
fim and 160 fim PACS channels and the 250 fim, 350 /mi, and 
500 fim SPIRE channels, over roughly 6 square degrees. The 70 
yum observations of Bl-E, however, are less sensitive due to the 
fast scan rate (60"/s) and low emission from cold material (see 
Section [3~TT i. Thus, we will not include the 70 fim observations 
in our discussion. For a full explanation of the observations, see 
Pezzuto et al. (2012, in preparation). 

The PACS and SPIRE raw data were reduced using HIPE 
version 5.0 and reduction scripts written by M. Sauvage (PACS) 
and P. Panuzzo (SPIRE) that were modified from the standard 
pipeline. For SPIRE reduction, we used updated calibration in- 
formation (version 4). The final maps were created using the 
scanamorphos routine developed by Roussel (2011, submitted 
to A&A). For consistency with other published maps, we set the 
pixel scal^H to the default size from the HIPE mapmaking tools, 
which results in pixels of sizes 6.4", 6.0", 10.0", and 14.0" 
for 160 fim, 250 fim, 350 fim, and 500 fim, respectively. We 



By default, scanamorphos sets a much smaller pixel scale. 



also adopt beam sizes of 13.4", 18.1", 25.2", and 36.6" for the 
160 urn, 250 yum, 350 ;um, and 500 um bands respectively (see 
Gri ffin et al J 120101: [Poglitsch et al . 1120101). Assuming a distance 
to Perseus of 235 pc dHirota et al.l 120081) . Herschel can detect 
structure on scales of ~ 0.04 pc at 500 fim. Figure Q] shows a 
three-colour image of Western Perseus with labels for Bl-E and 
other prominent subregions. 

Unlike previous ground-based submillimetre instru- 
ments, Herschel can detec t large-scale diffuse emission 
(e.g.. ISchneider et all 120101: iMiville-Deschenes et all 120101: 
lArzoumanian et all 1201 lb . Additionally, Herschel has excellent 
sensitivity to low-level flux. For example, Figure [2] compares 
SCUBA 8 50 urn observations of Bl-E (smoothed to 23" 
resolution; iDi Francesco et al.ll2008h with the new SPIRE 250 
fim observations (at 18" resolution). The SPIRE 250 fim data 
show prominent substructures not identified in the SCUBA 850 
fim data, though several faint features at 850 fim appear to agree 
with the brighter structures in the 250 fim map. With the higher 
sensitivity of Herschel, we are able to identify clearly structures 
in Bl-E that were too faint, i.e., < 3 <x, to be robust detections 
with SCUBA. 

2.2. GBT Observations 

We selected the nine highest peaks in our Herschel-denved H2 
column density map of Bl-E (see Section [372b for complemen- 
tary follow-up observations with the new K-band Focal Plane 
Array (KFPA) receiver on the GBT0 Our targets, hereafter 
named Bl-El to B1-E9 according to decreasing peak column 
density, were observed on 03 March 201 1 with single-pointings. 
We used the KFPA receiver with one beam and four spectral win- 
dows to observe each target simultaneously in NH3 (1,1), NH3 
(2,2), CCS (2i - 1 ), and HC 5 N (9-8) line emission at 23.6945 
GHz, 23.7226 GHz, 22.3440 GHz and 2 3.9639 GHz, respec- 
tively. Similar to lRosolowskv et al] (120081) . we made frequency- 
switching observations with 9-level sampling and 12.5 MHz 
bandwidth over 4096 spectral channels in each window to 
achieve a high velocity resolution of v c h ~ 0.0386 km s _1 at 
23.69 GHz. Bl-El to B1-E8 were observed for ~ 1 120 seconds, 
each. B1-E9 was observed for ~ 840 seconds. 

The GBT data were reduced using standard procedures in 
GBTIDI_Qfor frequency-switched data. In brief, individual scans 
from each spectral window were filtered for spikes or baseline 
wiggles, folded, and then averaged. Baselines were obtained for 
the averaged spectra by fitting 5th-order polynomials to line-free 
channel ranges at the low- and high-frequency edges of each 
band, and were then subtracted. To improve the detection levels, 
the data were smoothed with a boxcar kernel equal to two chan- 
nels in width. The reduced data were exported from GBTIDL 
into standard FITS files using routines of AIPS ++ developed by 
G. Langston. Further analysis was conducted in MIRIAE0 and 
IDL (Interactive Data Language). The 1 cr noise levels were typ- 



2 Due to the fast scan rate, the 160 fim beam is slightly elongated 
along the scan direction. Thus, we adopt the geometric mean (i.e., r = 
yfab) as the beam size , where the elongated beam dimensions are from 
iPoglitsch etaHd2010h . 

3 The GBT is an 100 m telescope operated by the National Radio 
Astronomy Observatory. 

4 GBTIDL is a special IDL package specific to the GBT. 

5 Multichannel Image Reconstruction, Interactive Analysis, and 
Display (MIRIAD) software is d eveloped by the Berkeley Illinois 
Maryland Array (BIMA) group. S ee l Saul t et alJd 19951) for more details. 
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Fig. 1. Three-colour image of the western half of the Perseus molecular cloud. Colour mosaic was generated using Herschel 160 
/mi, 250 yum, and 350 fim observations. The white box denotes our boundary for Bl-E for subsequent figures. The other prominent 
clumps in Perseus are also labeled. 



S B50 (Jy beam 1 ) S 250 (Jy beam l ) 

-0.1 0.1 0.2 0.5 1 1.5 2 




3 h 37 m s 40 s 20 s 36 m s 40 s 3 h 37 m s 40 s 20 s 36 m s 40 s 

a (2000) a (2000) 

Fig. 2. Observations of Perseus Bl-E fr om (a) SCUBA 850 /um a nd (b) SPIRE 250 fjxa maps. The SCUBA data came from the 
Extended SCUBA Legacy Catalogue (s eelDi Francesco et al.l20 08). Contours represent extinction levels of Ay — 5, 7, 8 magnitudes 
from the COMPLETE extinction map (Rid ge et al.l l2006a) and the filled circles show the beam sizes of 23" for smoothed SCUBA 
data at 850 fim and 18" for SPIRE data at 250 yum. 
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ically ~ 0.01 - 0.02 K. At 23.69 GHz, the GBT beam is ~ 33" 
and the beam efficiency is 77/, = 0.825. 

Table[T]gives the positions, peak H2 column density, and the 
detected line emission of these targets. NH3 (1,1) line emission 
was detected (> 3 <x) towards all nine targets and CCS (2; - lo) 
line emission towards several. The NH3 (2,2) and HC5N (9-8) 
lines were only detected toward B1-E2. 



3. Results 

3.1. SED Fitting to Herschel Data 

We corrected the arbitrary zero-poi nt flux offset in e ach H erschel 
band using the method proposed in Be rnard et alj d2010h that is 
based on a comparison with the P lanck HFI (DR2 version, see 
iPlanck HFI Core Team etai]|201 lb and IRAS data. In addition, 
each map was convolved to the resolution of the 500 fim map 
(36.6") and regridded to 14" pixels. The map intensities of the 
160 - 500 fim bands were then fit by the modified black body 
function, 

I v = k v B v {T)Y (1) 

where k y is the dust opacity, B v is the black body function, T is 
the dust temperature, and £ is the gas mass column density. Note 
that £ = /i»i//N(H2), where [i is the mean molecular weight, nin 
is the hydrogen mass, and N(H2) is the gas column density. For 
consistency with other papers in the Gould Belt Survey (e.g., 
lAndreet al.ll2010l) . 

A-^O.Kv/lOOOGHz/cir^g- 1 , (2) 

where B is the dust emissivity index. The SED fits were made 
using the IDL program mpfitfun by C. B. Markwardt. In brief, 
mpfitfun performs a least-squares comparison between the data 
and a model function by adjusting the desired parameters until a 
best fit is achieved. 

Fitting B requires a plethora of data, particularly along 
the Rayleigh-Jeans tail (e.g., A > 300 at 10 K), to re- 
move degeneraci es in the model fits (iDotv & Leun 2 [1992 
IShettv et al.ll2009h . Since we have only two photometric bands 
(350 //m and 500 //m) along the Rayleigh-Jeans tail, we as- 
sume B — 2, consistent with the adopted ft in other Herschel 
first-l ook studies (e.g., lAndre et all 120101 : lArzoumanian et all 
l201lb . The re is some evidence for B x 2 in recent Planck 
studies (see IPlanck Collaboration et al.ll201 1 alibi and references 
therein). Adopting B — 2 here provides good fits to our data (see 
below). 

Figure [3] shows the dust temperatures and column densities 
across Bl-E resulting from the modified black body fits to each 
pixel for the 160 - 500 //m bands, assuming flux uncertainties of 
15% based on calib ration uncertainties (see iGriffin et al.| [2010: 
Pogli tsch et al.ll20loh and B — 2. Both temperature and column 
density in Bl-E are highly structured (see Figures [3^ and[3t>), 
where regions of higher column density are also associated with 
slightly cooler temperatures and regions of lower column den- 
sity have slightly warmer temperatures. Most of the cooler, high 
column density material (> 6 x 10 21 cm -2 ) is clumped into a 
central ring-like structure (i.e., slight depression is seen towards 
the very centre). 

Figures [3J: and[3ji show the number histograms of tempera- 
ture and H2 column density, respectively. These distributions are 
non-Gaussian. The sample mean and standard deviation about 
the mean for the temperature and column density are 14.1 K±0.8 



K and (6.3 ± 2.7) x 10 21 cm 2 , respectively, where the standard 
deviation is computed from, 

where x is the population mean for the sample of size N. 

These mean values agree well with previously esti mated 
quantities for this region. For example, Sch nee et al.l ^2005) mea- 
sured a mean dust temperature of ~ 14 K for the Bl-E region 
using the ratio of IRAS 60 /mi and 100 fim flux densities to es- 
timate dust temperature. Additionally, extinction data from the 
COMPLETE survey (see contours in Figures [2] and 01 suggest 
a mean column density of (5.3 ± 1.5) X 10 21 cm -2 , assuming 
N(H2)/Ay = 10 21 cm -2 mag -1 . Using a higher resolution ex- 
tinction map from S. Bontemps, we find a mean column density 
of (6.3 ± 3.6) x 10 21 cm -2 . These two column densities agree 
very well with our measured value of (6.3 + 2.7) x 10 21 cm" 2 
and they are also similar to the threshold column density for 
dense core form ation from the extinction analysis in Perseus by 
iKirket all (120061) . 5 x 10 21 cm 2 . For comparison, Figure |5Jl in- 
cludes the Kirk et al. column density threshold as a dashed line. 
Although Bl-E does not appear filamentary, there is ample ma- 
terial frojn_which densecores may form. With higher resolution 
data, Andre et"a"D d201 ll) found a threshold column density of 
7 X 10 21 cm" 2 for dense structures to form in Aquila via thermal 
instabilities along a filament at 10 K. 

3.2. Column Density Profiles 

One of the more prominent core models is the Bonnor-Ebert 
sphere dBonnorll 1956b lEbertlll955l) . which represents the den- 
sity profile of a sphere in hydrostatic equilibrium under the in- 
fluence of an external pressure. The Bonnor-Ebert sphere has 
a flat inner density distribution and a power-law density down- 
trend at larger radii. Many prestellar cores, i.e., dense cores that 
are gravitationally bound but do not show evidence of a central 
luminous protostar, have shown Bonnor-Ebert- like profiles (e.g., 
I Ward-Thompson et al.lll999UAlves et~aDl200lh . 

With the excellent spatial resolution of Herschel, we can 
measure the column density profiles for individual Bl-E sub- 
structures. First, we identified the locatio ns of peak column den - 
sity using the 2D Clumpfind algorithm (Willia ms et al.l fl994>. 
Briefly, Clumpfind identifies intensity peaks and then uses 
closed contours at lower intensity levels to assign boundaries. 
We will discuss the nine highest column density substruc- 
tures identified by Clumpfind (Bl-El to B1-E9) for this paper. 
Second, we measured the azimuthally-averaged column density 
profile of our nine sources using the ellint task in MIRIAD. For 
simplicity, we used circular annuli of 10" width for r > 1" . For 
the central radii (r < 7"), we assume the peak column density. 
We caution that several of our sources appear elliptical and our 
circular approximation is meant to provide a broad, first-look 
analysis. 

Figure [4] shows the column density profiles, where the col- 
umn density values are plotted from the centre of each annu- 
lus. For comparison, Figure|4]also shows the beam profile (solid 
grey curve), a "generic" column density profile convolved with 
the beam (dashed grey curve) and an analytical Bonnor-Ebert 
profile convolved with the beam (dotted curve). For both ana- 
lytic p rofiles, we follow the approximation from Dapp & Basu 
(2009) and assume a temperature of T — 10 K, a central density 
of n = 10 6 cm" 3 and a constant of proportionality of k = 0.54, 
where k at 0.4 for a singular isothermal sphere and k as 1 
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Table 1. GBT Target Information 



Source 


a (J2000) 


6 (J2000) 


Peak N(H 2 ) fl 


Detections 6 




(h:m:s) 


(° ' ") 


(10 21 cm- 2 ) 




Bl-El 


3:35:55.0 


31 


14:16 


18.5 


NH 3 (1,1), CCS (2i-l ) 


B1-E2 


3:36:04.4 


31 


11:47 


16.7 


NH 3 (1,1), NH 3 (2,2), CCS (2, - 1 ), HC 5 N (9-8) 


B1-E3 


3:35:52.1 


31 


15:39 


15.1 


NH 3 (1,1), CCS (2, - 1 ) 


B1-E4 


3:35:51.0 


31 


12:34 


14.9 


NH 3 (1,1) 


B1-E5 


3:36:37.3 


31 


11:41 


14.5 


NH 3 (1,1) 


B1-E6 


3:36:41.0 


31 


15:05 


13.7 


NH 3 (1,1) 


B1-E7 


3:36:39.0 


31 


14:27 


13.6 


NH 3 (1,1), CCS (2, - 1 ) 


B1-E8 


3:36:05.0 


31 


14:28 


13.1 


NH 3 (1,1), CCS (2,-l ) 


B1-E9 


3:36:18.5 


31 


14:31 


12.9 


NH 3 (1,1) 



Notes. (fl) Peak /feri'c/ieZ-derived H 2 column density towards the sources. See Section[3]for more details. <6) Line emission detected towards each 
source. 




Fig. 3. SED-fitting results for Bl-E. Top panels shows the temperature (a) and H2 column density (b) maps across Bl-E as measured 
from SED-fitting to 160 - 500 yum data, according to Equations Q] and [2] and assuming B = 2. Bottom panels show the histograms of 
the above maps. The temperature histogram (c) uses a bin size of 0.3 K and the H2 column density histogram (d) uses a bin size of 
1 x 10 21 cm 2 . For comparison , the d ashed line in panel (d) indicates the observed column density threshold (from extinction) for 
core formation from Kir k et all (120061) . 



for a collapsing cloud. Our observed column density profiles 
are much wider than the models, likely due to contaminating 
material in the foreground or far-background along the line of 
sight, hereafter called LOS material. Towards the centre of each 
source, the analytic models and observed profiles are more sim- 
ilar, since the source column density dominates over the LOS 



level, whereas in the power-law roll-off, the LOS material is 
likely more significant and the profiles deviate from the analytic 
models. Unfortunately, the column density of such extended ma- 
terial is difficult to distinguish from the source column density. 
To estimate its contribution towards each source, we used the 
subsequent column density at the location where the power-law 
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slope in the column density profile flattened or began to increase. 
Thus, the LOS material level ranges from ~ 8 x 10 21 cm -2 for 
B1-E9 to ~ 11 x 10 21 cirT 2 for Bl-El or roughly 60% of the 
peak column density. These LOS column densities are very con- 
servative and could overestimate their actual contributions by as 
much as a factor of 2. Indeed, these values are larger than the 
mean Bl-E column density of 6.3 x 10 21 cirr 2 . 

Figure [5] shows a comparison of our nine column density 
profiles after correcting each profile for LOS material and nor- 
malizing to the respective peak column density, and includes 
the Bonnor-Ebert profile, generic profile, and beam profile from 
Figure |4] All nine profiles follow a shape similar to those of the 
analytic models, with a flat centre and steep falloff towards larger 
radii. Additionally, all profiles but B1-E5 generally appear more 
centrally concentrated than the models, implying that B1-E5 is 
the least compact. The column density profiles do show some 
differences with the models, however, such as steeper fall-offs. 

3.3. Substructures 

We estimated source sizes and masses based on Gaussian fits 
to the LOS-subtracted column density profiles. Source size was 
defined by the FWHM of the Gaussian (i.e., R = FWHM/2). 
Since B1-E6 and B1-E7 have a small projected separation of 
~ 50", we truncated both to 31" radiQto ensure that we mostly 
measured column densities unique to each. Thus, our analyses 
of B1-E6 and B1-E7 are biased toward the denser, central re- 
gions unlike those of the other sources. Subtracting out the LOS 
material limits the influence from nearby extended material and 
generally fits only the denser, embedded object. If we include the 
LOS material, the estimated source sizes generally increase by 
< 20%. Isolated substructures, like B1-E2 or B1-E8, have sizes 
that vary by ~ 5% and substructures embedded in extended pro- 
files like, Bl-El or B1-E4, have sizes that vary by 24% and 37%, 
respectively. 

Table [2] lists the adopted properties for our sources, includ- 
ing their deconvolved radii, estimated (upper and lower) masses, 
and average densities, assuming they are perfect spheres and di- 
viding the (upper and lower) source masses with their spherical 
volumes, f tt^, where Rj is the deconvolved radius. We con- 
sider two mass and density limits: the lower limits subtract out 
the LOS material and the upper limits include the LOS mate- 
rial. Note that, by our definition of source size, we are measur- 
ing the inner regions of each profile where core precursors are 
most likely to arise. Since we are using the Gaussian FWHM 
to estimate the source size, the true source sizes, including any 
diffuse envelopes, may be larger by a factor of 2. Source mass 
was estimated by summing over the column density and assum- 
ing a mean molecular weight of fi = 2.33. These results provide 
a reasonable first look at the relative sizes and masses of these 
sources. 



3.4. Line Emission 



Table 2. Column Density Determined Properties 



iKirk et al.l (120071) and iRosolowskv et all d2008l) each attempted 
to identify dense gas towards Bl-E via high spatial resolution 
line observations in N2H + and NH3, respectively. They observed 
several "blind" pointings towards Bl-E and found no strong de- 
tections. Unlike previous continuum studies, our Herschel data 
have identified several cold, dense substructures in Bl-E (see 
Figure [3]) and these were not directly probed by either Kirk 
et al. or Rosolowsky et al. Using our Herschel results, we se- 



Source 


R/' 


M" 


rf 






(AU) 


(Mq) 


(10 4 cm" 3 ) 


Bl-El 


7.1 xlO 3 


0.6 - 1.6 


5.9- 


17 


B1-E2 


6.9 xlO 3 


0.6 - 1.4 


6.7- 


15 


B1-E3 


7.7 xlO 3 


0.5 - 1.5 


4.0- 


12 


B1-E4 


6.0 xlO 3 


0.3- 1.0 


4.3- 


18 


B1-E5 


9.3 xlO 3 


0.5 - 2.0 


2.4- 


9.1 


Bl-E6 rf 


5.9 xlO 3 


0.3 - 1.0 


5.7- 


17 


Bl-E7 rf 


5.9 xlO 3 


0.3-0.9 


5.6- 


17 


B1-E8 


7.2 xlO 3 


0.3 - 1.2 


3.1 - 


12 


B1-E9 


5.4 xlO 3 


0.3-0.8 


5.8- 


18 



Notes. <n) Object radii were estimated from Gaussian fits to the col- 
umn density profile (R = FWHM/2). Radii were deconvolved with a 
36.6" beam. (A) Masses are taken from integrating out the column den- 
sity profiles. The lower limit measurements subtract the LOS material 
from the column density profiles. The upper limit measurements do not 
subtract the LOS material from the column density profiles. {c) Average 
density determined from the mass and volume. <<fl B 1-E6 and B 1-E7 are 
truncated to measure column density unique to each source. 



lected the nine substructures with the highest column densities 
for follow-up observations with the GBT in NH3 (1,1), NH3 
(2,2), CCS (2j - 1 ), and HC 5 N (9-8) line emission. Figure 
shows the locations of our nine targets (see also TableQ}. 




31 06 



3 n 37 m s 



■15 



30 s 15 s 36 m s 
Q (2000) 



45' 



30' 



Fig. 6. Locations of substructures in Bl-E overlaid on our 
//eric/ie/-derived column density map. Numbers indicate the po- 
sitions of peak column densities and the relative magnitude of 
the column density peak (see also Table Q]). The thin black cir- 
cles show the deprojected source sizes (see Section I3.3I ). The 
dark filled circle shows the GBT beam (~ 33") at 23.69 GHz. 



3.4.1. NH 3 

We detected NH3 (1,1) emission towards all nine targets. 
We fit the NH3 (1,1) hyperfine lines with multiple Gaussian 
components using the following equation; 



Deconvolved radii are ss 25" 



18 

t(v) = r u ^a ; exp 

1=1 



-4 In 2 



I v-x'i - VfaA 2 
I Av J 



(4) 
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Fig. 4. The azimuthally-averaged column density profiles for the nine sources in Bl-E, assuming a distance of 235 pc. The profiles 
were measured using circular annuli with a thickness of 10" for r > 7". Note that the central area is defined by a circle of radius 
7" but plotted at 3.5". The dashed line indicates our estimate for the non-source line-of-sight (LOS) material. For comparison, the 
dotted curve and dashed curve illustr ate a Bonnor-Ebert column density profile and a "generic" [1 + (r/a) 2 ]~°- 5 column density 
profile followin dDapp & Ba su (2009). For both analytic curves, we assume a temperature of 10 K and a central density of n — 10 6 
cm -3 and convolved the profiles with a 36.6" beam. The grey solid curve shows the beam profile. The analytical profiles and the 
beam profile are scaled to the peak column density. Note that both axes are logarithmic. 



where T14 is the optical depth for the (1,1) transition, a, is the 
transition weight, v, is the velocity of the hyperfine line (for 
a given rest frequency), vi sr is the velocity of the source with 
respect to the local standard of rest, and Av is the FWHM of 
the line. We use d th e hyperfine frequenci es and weights from 
iKukolichl d 19671) and Rvdbe ck et all (1 19771) and derived the val- 
ues for T11, v^, and Av from the fits. Figure [7] shows example 
fits to the NH3 spectra. Table [3] lists centroid velocity (v/ sr ) and 
uncorrected velocity line width (Av) from the NH3 (1,1) fits, the 
kinetic gas temperature from NH3 (1,1) and (2,2) line compar- 
isons, and the resulting thermal and non-thermal velocity disper- 
sions for all nine sources. 

For B1-E2, we could derive the kinetic gas temperature us- 
ing NH3 (1,1) and NH3 (2,2) line emission, whereas we could 
only determine upper limits for our other sources. We use our 
NH 3 -derived kinetic temperature of T K = 10.27 K + 0.38 K 
for B1-E2 and adopt the mean NH3 -derived kinetic tem pera- 
ture of 1 1 K from the survey of iRosolowskv et alJ (l2008h for 



our other sources. As an estimate of the uncertainty, we use the 
dispersion in th e kinetic temperatu res from the Rosolowsky et al. 
data, 2 K, (see En och et 101 12008). Note that these gas tempera- 
tures correspond to the denser, colder inner layers where ammo- 
nia is excited, whereas the dust temperatures (see Section 13. U 
represent an average of all the material along the line of sight, 
i.e., including the lower density envelopes. For more informa- 
tion lon^xttacting kinetic temperatures from NH3 line emission, 
see iFriesen et aild2009l) . 

Table[3]also lists the thermal velocity dispersion, cr T , and the 
non-thermal velocity dispersion, cr NT . We determined crj from 
the kinetic temperature and cr NT from the line widths via, 



0~NT 



hT K 



y f*NH 3 



m H 



■ 0~T 
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Fig. 5. Normalized column density profiles for sources in Bl-E, assuming a distance of 235 pc. Column density was determined 
from azimuthally-averaged annuli (see Figure[4]i. Furthermore, the LOS material was estimated for each source and subtracted from 
the azimuthal average. The grey solid, dotted, and dashed curves are the same as for Figure|4] 

Table 3. Properties from NH 3 (1,1) and NH 3 (2,2) Line Spectra 



Source 


Vhr 


Av 


1 


T K " 






<r nt 




(km s" 1 ) 


(kms- 1 ) 


(K) 


(kms- 1 ) 


(km s" 1 ) 


Bl-El 


7.61 ± 0.02 


0.994 ± 


0.055 


< 27 


0.073 ± 


0.007 


0.415 ± 0.025 


B1-E2 


7.67 ± 0.002 


0.293 ± 


0.005 


10.27 ± 0.38 


0.071 ± 


0.001 


0.101 ± 0.003 


B1-E3 


7.47 ± 0.01 


0.681 ± 


0.030 


< 16 


0.073 ± 


0.007 


0.280 ± 0.015 


B1-E4 


7.82 ± 0.01 


1.080 ± 


0.039 


< 27 


0.073 ± 


0.007 


0.452 ± 0.018 


B1-E5 


7.57 ± 0.01 


0.832 ± 


0.022 


< 19 


0.073 ± 


0.007 


0.344 ±0.011 


B1-E6 


7.61 ±0.01 


0.905 ± 


0.032 


<23 


0.073 ± 


0.007 


0.376 ± 0.015 


B1-E7 


7.31 ±0.01 


0.708 ± 


0.032 


< 24 


0.073 ± 


0.007 


0.290 ± 0.016 


B1-E8 


6.78 ± 0.01 


0.537 ± 


0.021 


< 20 


0.073 ± 


0.007 


0.213 ±0.012 


B1-E9 


7.11 ± 0.01 


0.688 ± 


0.038 


< 32 


0.073 ± 


0.007 


0.281 ±0.018 



Notes. (a) The velocity FWHM from the best-fit models uncorrected for channel resolution. {b) Kinetic temperature towards each source. For B 1 -E2, 
we measured the kinetic temperature from NH3 (1,1) and NH3 (2,2) line emission. For our other sources, we used a 3 a upper limit for NH3 (2,2) 
to estimate the upper limit of T K . (c) The thermal line width of B 1-E2 was determined using the d erived kinetic temperatur e, T K = 10.27 K ± 0.38 
K. For the remaining sources, we adopted 1 1 K, the mean NH3 -deriv ed kinetic tempera ture from Rosolows kv et alJ {2008) with an uncertainty of 
2 K from the dispersion of the kinetic temperature in that survey (see Enoch et al. 200j|). 



where kb is the Boltzmann constant, Tk is the kinetic temper- 
ature, jUwh, = 17.03 is the mean molecular weight of NH3, 
m.H is the atomic hydrogen mass, and <T b s is the total cor- 
rected velocity dispersion of the lines, cr b s — Av r / V81n2 for 
Av c 2 = Av 2 - v 2 ;i . Using the NH3 study by iRosolowskv et alJ 
(2008J), we adopted T K = 1 1 K for all our targets except B1-E2, 
where we used our derived gas temperature, T K = 10.27 K+0.38 
K. 



3.4.2. CCS and HC 5 N 

CCS (2i - lo) was detected towards five targets (see Table [TJ 
and HC 5 N (9-8) was detected towards only B1-E2. Figure [8] 
shows examples of these spectra. We fit the CCS {2\ - lo) and 
HC5N (9-8) spectra with single Gaussians. The velocities deter- 
mined for most CCS {2\ - lo) detections and the single HC5N 
(9-8) detection were similar to those found from NH3, suggest- 
ing these dense gas tracers are probing the same region as the 
NH3 emission. The CCS (2\ - lo) spectra for Bl-El and Bl- 
E8, however, revealed fairly different values. For Bl-El, the 
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Fig. 7. Fits to NH3 (1,1) and (2,2) spectra. Note the hyperfine components visible for the NH3 (1,1) spectrum of B1-E2. The NH3 
(2,2) spectra were smoothed by a boxcar with a 2-channel width to improve the visual appearance. For comparison, we show two 
undetected NH3 (2,2) observations. 



CCS (2i - lo) line emission appears redshifted by ~ 0.3 km s 
(14 cr) with respect to the NH3 (1,1) line emission, whereas for 
B1-E8, the CCS (2i - lo) line emission is redshifted by ~ 0.1 
km s _1 (1.5 cr). The redshifted CCS (2\ - lo) emission in Bl- 
El may be tracing a different region of dense gas fro m the NH3 
emiss ion, such as a region undergoing infall (e.g., see Swi ft et all 
I2005h . 3 

4. Discussion 

Our Herschel data have revealed extensive substructures within 
Bl-E for the first time. Such faint substructures may represent 
an early evolutionary stage in star formation, where a first gen- 
eration of cores is forming from lower density material. Thus, 
it is important to characterize the substructures seen in these 
Herschel maps. 



4.1. Comparison with Jeans Instability 

The minimum length scale for gravitational fragmentation in a 
purely thermal clump or cloud is described by the Jeans length, 
Aj - (c 2 s ji IGp) 1 ! 2 , where c s is the thermal sound speed, G 
is the gravitational co nstant, and p is the density of material 
(IStahler & Pallal2005l) . In terms of temperature and number den- 
sity, we find, 

/ T \ - 5 / n r - 5 

^ - 21 (K)l) (low) pc - (?) 

Ideally, the Jeans length should be calculated for the gas temper- 
ature and gas density. In small, dense scales (i.e., > 10 s cm -3 ; 
iDi Francesco et al.l20 07). the gas and dust temperatures are cou- 
pled, whereas in larger, more diffuse scales, the gas a nd dust are 
decoupled and the gas may be warmer than the dust (Young et al. 



9 



Sadavoy et al.: Herschel Observations of Perseus Bl-E: A First Look at Core Formation 



0.10 
0.05 
0.00 

-0.05 

0.1 
0.08 
0.06 
0.04 
0.02 
0.00 
-0.02 
-0.04 



CCS (2-1) - 



1 hill I 



5 10 
v (km/s) 



1 5 20 



- B1-E3 

In! Ill 


llllll 




CCS (2- 
Li ||, 


1) - 
I I lit 

lil ? 


Ml ft 


rijj 




f 


I I' > 
|'1 


'I'll 

f H 




0.6 r 

0.5 I 
0.4 j- 
0.3 : r 
0.2 r 

0.1 i- 

0.0 



CCS (2-1) 




i will ' 



5 10 
v (km/s) 

Fig. 8. Fits to the detected CCS (2i - 1 ) and HC 5 N (9-8) spectra. The weaker CCS (2i - 1 ) spectra and the B1-E2 HC 5 N (9- 
8) spectra were all smoothed by a boxcar with a 2-channel width to improve the visual appearance of the lines. The B1-E2 
CCS (2i - lo) spectrum was not smoothed. 



12004 Pi Francesco et a D2007t lCeccarelli et al. 2007). The cen- 
tral region of Bl-E containing the substructures, i.e., the area 
contained within N(H2) > 6 x 10 21 crrT 2 , has a mass of ~ 110 
M and an effective radius (R e ff = y/A/n) of ~ 0.46 pc. 
Approximating Bl-E as a sphere, we find an average density 
of ~ 5 x 10 3 cm 4 , suggesting that the the gas and dust in Bl- 
E are likely decoupled such that the mean dust temperature of 
~ 14 K (see Section [37Tb is a lower limit to the gas temperature. 
Assuming a lower limit gas temperature of 14 K and n — 5 x 10 3 
crrT 3 , the lower limit Aj * 0.35 pc. Note that the kinetic tem- 
peratures measured in Section 13.4.11 correspond to the smaller, 
denser scales where the gas and dust temperatures should be cou- 
pled and not the scales measured here. 

The minimum projected separations between our nine sub- 
structures are generally between 0.1 pc and 0.19 pc, with a me- 
dian value of ~ 0.13 pc for all nine sources. Recall, however, 
that B1-E6 and B1-E7 have an angular separation of ~ 0.05 pc, 
a factor of 2 closer than any other pair of substructures. Despite 
their close proximity, B1-E6 and B1-E7 show different centroid 
velocities (see Table |3), and therefore, cannot be considered a 
single object. Also, we can only measure the projected 2-D sep- 
arations between our nine substructures and not their physical 



3-D separations. Assuming a typical inclination angle of 60°, 
our median separation is ~ 0.15 pc for all nine substructures or 
~ 0.17 pc excluding B1-E6 and B1-E7. 

The median minimum separation is a factor of two smaller 
than our best estimate of the Jeans length, though our Jeans 
length estimate is, itself, uncertain within a factor of a few. 
The fact that we find a minimum separation that is less than 
the Jeans length is still interesting and could indicate gravita- 
tional fragmentation followed by bulk contraction of the Bl- 
E group, but the overall uncertainties in determining the mini- 
mum separations and Jeans length make this difficult to validate. 
Additionally, Bl-E does not appear to be filamentary (see Figure 
[6]l and therefore, the clump could be extended along the line of 
sight resulting in more significant radial distances between sub- 
structures. 

For the substructures, we use the Jeans mass to explore 
their thermal stabilities. Similar to the Jeans length, the Jeans 
mass is a critical m ass scale for gravitational fragmentation 
(IStahler & Pallall2005h . This critical mass can be described as, 
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The Bl-E substructures have average densities of n ~ 1 x 10 s 
crrT 3 (see Section l3~3l >. Furthermore, we found a gas tempera- 
ture for B1-E2 of ~ 10 K. Note, that this gas temperature cor- 
responds to the denser B1-E2 substructure and should not be 
assumed for the lower density Bl-E clump. Assuming T = 10 
K and n = 1 x 10 5 cirT 3 , the critical Jeans mass is My * 0.9 
M Q . Several of the substructures have mass limits on order of 
the Jeans mass (see Table suggesting that these sources are 
approaching a critical, unstable mass. Thus, the Bl-E substruc- 
tures are interesting prospects for future studies of substructures 
evolution. 



4.2. Time Scale for Interactions 

From the Herschel observations, the Bl-E substructures appear 
linked (see Figure [6), which might indicate that they are inter- 
acting. The timescale for any two substructures to interact can 
be estimated from the typical separation between the substruc- 
tures and their velocity dispersion. For simplicity, we adopt the 
average projected separation between all nine substructures of 
~ 0.40 pc. For an inclination angle of 60°, the typical separa- 
tion between the substructures is d ave ~ 0.46 pc. To estimate 
the velocity dispersion of the sources, we use the dispersion of 
the centroid velocities of NH3 (1,1) from all nine objects. The 
weighted ID centroid velocity dispersion is ct\d ~ 0.24 km s _I 
and assuming symmetry, the 3D velocity dispersiorQ is cr 3D ss 
0.42 kms~'. Thus, the Bl-E substructures have an interaction 
timescale of ~ 1 Myr, a factor of two larger than recent estimated 
prestellar core lif etimes and less than the expected protostellar 
core lifetime (e.g iKirk et al.l2005l:IWard-Thompson et alj2007t 
lEnoch et"aT]|2008l) . Furthermore, this timescale corresponds to 
the interaction time between any two substructures and not 
the interaction time between nearest neighbour substructures. 
Given that the nearest neighbour separation is a factor of ~ 3 
smaller than the typical separation adopted here, these substruc- 
tures have the potential to interact prior to forming stars. Thus, 
competitive accreti on may be significant during further evolu- 
tion of Bl-E (e.g.. iBonnell et all l200lt iKramholz etail [2005 ; 



Table 4. Virial Analysis Properties 



IBonnell & Bate 2006). 



4.3. Comparison with Virial Equlibrium 

Not all dense substructures will form cores and stars. For ex- 
ample, a core with kinetic energy from internal motions that ex- 
ceeds its gravitational binding energy may be transient and un- 
able to form p ersistent dense structures or stars. According to 
iMcKed (|T999), a core is gravitationally bound if the virial pa- 
rameter, or, is 



M virial 5o-t,R 
a = — — — = < 2 



M 



GM 



(9) 



where cr v is the velocity dispersion, R is the radius, and M is the 
mass. 

Bl-E has an effective radius of 0.46 pc (see Section PTTl i. 
We estimate the velocity dispersion using the centroid velocity 
dispersion from Section H~2l <x 3D = 0.42 km s _1 . This velocity 
dispersion was derived from the motions of the individual sub- 
structures and thus, reflects the turbulence in the cloud that first 
created the substructures rather than the thermal velocity dis- 
persion of the gas. The expected thermal velocity of the gas is 
less than 0.42 km s _1 , and thus non-thermal support is impor- 
tant. Assuming R = 0.46 pc and cr v = 0.42 km s _1 , Bl-E has a 



Source 


(km s" 1 ) 


£>"„/' 


„ a 

"S 


Mvirial 

(M ) 


a 




Bl-El 


0.460 ± 


0.030 


2.3 + 


.0.3 


8.5 ±4.8 


5- 


14 


B1-E2 


0.216 ± 


0.005 


1.1 ± 


0.03 


1.8 + 0.9 


1 - 


■3 


B1-E3 


0.342 ± 


0.022 


1.7 ± 


: 0.2 


5.1 ±2.9 


3- 


10 


B1-E4 


0.493 ± 


0.024 


2.5 ± 


0.3 


8.2 ±4.5 


8- 


27 


B1-E5 


0.397 ± 


0.018 


2.0 + 


:0.2 


8.3 ±4.5 


4- 


17 


B1-E6 


0.424 ± 


0.022 


2.2 + 


0.2 


6.0 ±3.3 


6- 


20 


B1-E7 


0.351 ± 


0.023 


1.8 + 


0.2 


4.1 ±2.3 


5- 


14 


B1-E8 


0.291 ± 


0.021 


1.5 + 


:0.2 


3.4 ± 2.0 


3- 


11 


B1-E9 


0.343 ± 


0.025 


1.7 + 


:0.2 


3.6 ±2.1 


4- 


12 



Assuming a\ D 



of + o\ + a* 



3<T 



Notes. (a) Velocity dispersion divided by sound speed, assuming T = 1 1 
K for all substructures except B1-E2 where T = 10.27 K. (w Virial pa- 
rameter considering the upper and lower mass limits for each substruc- 
ture (see Table 



virial mass of ~ 96 M , or a « 0.9 for a clump mass of 1 10 M Q 
(see Section l4~TT i. These results imply that Bl-E is itself gravita- 
tionally bound and thus, likely to form dense cores and stars in 
the future. 

For the individual substructures, velocity dispersions were 
determined from NH3 (1,1) line widths, after estimating the non- 
thermal component {(Jnt', see Table [3j> and the thermal compo- 
nent (c s ) for a mean molecular weight of /jh 2 = 2.33. Table |4] 
shows our values for the total velocity dispersion, <r v , the ra- 
tio of the velocity dispersion to the thermal sound speed, the 
virial mass,M,,, nfl ;, and the virial parameter, a, for each substruc- 
ture. Most of the virial parameters are » 2, suggesting that these 
substructures themselves are not gravitationally bound and thus, 
may not necessarily represent persistent objects. B1-E2 has the 
smallest virial parameter, 1 < a < 3, however, and may be it- 
self gravitationally bound. B1-E8 and B1-E3 have lower limit 
estimates of a ~ 3 and given that we know a only within a fac- 
tor of a few, these sources may also be gravitationally bound. 
These virial limits, however, depend greatly on the source mass, 
which can vary by a factor of ~ 4 depending on our treatment 
of the LOS material (see Section [33T l. To be bound, these candi- 
dates must have negligible LOS material towards them, which is 
very unlikely. Thus, B1-E2 is the only strong candidate for being 
gravitationally bound. 

Most Bl-E substructures have very large virial parameters 
and are thus, expected to be unbound. In comparison, the Jeans 
mass (see Section PTTT i is a mass scale for only gravitational frag- 
mentation in a thermal clump and does not include the contri- 
butions from non-thermal support. Several of the substructures 
with a » 2 have masses on order of the Jeans mass, ~ 0.9 M Q , 
suggesting that the Jeans mass analysis is too simplistic and does 
not well represent the evolutionary state of these structures. For 
example, the substructures may have significant turbulent mo- 
tions. With only one potentially bound substructure (B 1 -E2), it is 
difficult to determine if Bl-E will eventually form several cores 
or no cores at all (i.e., these substructures are all transient). 

The substructures themselves show large line widths (see 
Table atypical of dense prestellar cores, which are subsonic. 
In particular, B1-E2 has the narrowest line profiles (by at least a 
factor of 2) and the strongest line detections, suggesting that this 
object is the most evolved of the substructures. These observa- 
tions hint at a possible dynamic evolution, where dense cores 
first contain significant non-thermal, turbulent motions which 
dissipate into coherent, quiescent dense cores. If so, the sub- 
structures in Bl-E may become more bound over time as their 
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turbulent support is dissipated, a process which may explain the 
narrower lines toward B1-E2. 

Note that we assume Tk - UK for our substructures (ex- 
cept B1-E2). In Section 13.4.11 we found upper limits for kinetic 
temperature > 1 1 K (see Table [3}. Furthermore, we are using 
a simplified virial equation in Equation djj). Ideally, the virial 
equation includes magnetic pressure and surface pressure terms 
which can either assist gravit y in the collapse or h elp in the sup- 
port of a clump (e.g., see iNakamura & Ll 12008). Our observa- 
tions, however, do not measure these quantities and thus, we use 
a simplified virial equation. Nevertheless, B1-E2 is special with 
respect to the other substructures. Indeed, we may have the first 
observations of a core precursor. 



4.4. Comparison with Other Star Forming Regions 

Bl-E has column densities well abo ve ~ 5 X 10 21 cix r 2 , the 
core formation threshold reported by Kirk et all (120061) . Since 
the other clumps with similarly high column densities are ac- 
tively forming stars, it is reasonable t o expect B l -E to do the 
same in the near future. Furthermore, lLada etalJ (1201 Oh found 
a good correlation between the number of YSOs in a cloud (or 
clump) with the cloud mass above Ak > 0.8 (or Ay > 7.3) for 
several star forming regions. Based on their extinction threshold, 
we find a clump mass of ~ 90 M for N(H 2 ) > 7.3 x 10 21 crrT 2 
and thus, we could expect ~ 10 YSOs to form in Bl-E. With 
Herschel, we have detected several substructures with moderate 
densities of ~ 10 5 crrT 3 (see Tabled, s imilar to what is often d e- 
fined for dense cores, i.e., > 10 4 crrT 3 (Ber gin & Tafa lla 2007). 

These substructures, however, were not well observed by 
SCUBA at 850 pm or Bolocam at 1 mm (see Figure and 
thus, are not dense cores in the traditional definition. The 
substructures in Bl-E also have supersonic velocity disper- 
sions. In general, lo w-mass starless cores have subsonic veloc- 
ity di s persions (e.g iMverslI 19831: iKirk et al]l2"o"o"7l lAndre et all 



spersions (e.g ^ 
2007; iPineda et al.l 120101) whereas prot ostellar cores ca n have 
supersonic velocity disp ersions (e.g., iGregersen et al.l 119971; 
iDi Francesco et al1l200ll) . Since we see no evidence of proto- 
stellar activity whatsoever (i.e., from Spitzer data), these sub- 
structures are likely unevolved. 

Thus, the substructures in Bl-E are likely core precursors 
and more observations of this region should provide some in- 
sight into core formation (transient or persistent) in molecular 
clouds. In particular, Bl-E is isolated from YSOs and young 
stars, so the region is unaffected by environmental processes 
such as outflows or winds. A near by expanding dust s hell, seen 
in IRAS 60 pm and 100 fxm data dRidge et al.ll2006bl) . may be 
sweeping material towards Bl-E, though the ring itself does not 
appear to be directly interacting with this clump. This shell may 
have a larger impact on the evolution of structure in B 1 -E in the 
future. 

A relatively pristine core forming region is exceedingly 
rare. The only additional case may be t he inactive high extinc- 
tion clump (211) identified in L1495 bv lSchmalzl etalJ (l2010l) . 
though further detailed studies of that clump are still necessary to 
determine what, if any, substructures are forming cores and stars. 
A collection of starless cores in Aquila form a similar grouping 
as in B 1 -E, however this region has two kno wn protostars nearby 
(< 1 pc) and so is not as pristine as Bl-E dKonvves et al.l f2010, 
see their Figure 5). 



4.5. Comparison with Core Formation Models 

Core formation simulations attempt to recreate molecular 
cloud fragmentation, implementing factors such as gravita- 
tional collapse, magnetic pressure, and turbulent energy. Many 
recent studies have argued that models of core formation 
must consider both ambipolar diffusion and turbulence (e.g. , 
| Li & Nakamural 12004 iDib et alJ 120071: INakamura & Lil [2008; 
iBasu et al.ll2009allbl:lGong & Ostrikerll201 ll) . Magnetic fields in- 
troduce an additional pressure support that opposes gravitational 
collapse, whereas turbulence either opposes collapse by intro- 
ducing kinetic energy support or induces collapse via energy dis- 
sipation from compressive shocks. 

Several recent studies have examined the effects of dif- 
ferent magnetic field and turbulent flow properties (e .g., 
iLi & Nakamurall2004tlBasu et al.ll2009cJ:lPrice & Batell2009l) . In 
particular, these simulations examine fragmentation in subcrit- 
ical (where magnetic fields dominate over gravity) and super- 
critical (where gravity dominates over magnetic fields) regimes 
with dissipating turbulent flows. In the subcritical simulations, 
the fragmenting clouds generally form cores in relative isolation 
and with less filamentary morphologies , such as ring-like dis - 
tributions of loose core groups (see also ILi &~Na kamura 2002). 
These cores tend to form more slowly, i.e., fragmentation is sup- 
pressed by the strong magnetic field, and since large-scale tur- 
bulent flows are gen erally damped, the velocity field is generally 
subsonic. (Note that Basu et al. 2009a, find that strong magnetic 
fields can induce oscillations and create a supersonic velocity 
field, i.e., < 3c,.) In contrast, supercritical simulations gener- 
ally form cores associated with filaments, and these cores tend 
to form relatively quickly and can move supersonically through 
their environment. 

An isolated core forming region is necessary to test model 
predictions and constrain initial conditions of dense gas before 
the onset and influence of nearby stars and internal protostars. 
Thus, Bl-E is an ideal target for such measurements (see Section 
14.41 ). From our observations, Bl-E contains a small, loose group- 
ing of substructures. Although the substructure morphology ob- 
served in Bl-E is not entirely ring-shaped, the distribution hints 
at an inclined ring-like structure. Furthermore, if star formation 
is triggered, it is surprising that Bl-E, which is bookended by 
two highly active and more evolved clumps, IC 348 to the far 
east and NGC 1333 to the far west, contains no evidence for 
evolved star formation itself. Furthermore, IC 348 shows many 
evolved pre-main sequence star s and NGC 1333, though appear- 
ing younger, has many YSOs dBallv et all 120081; lHerbstl T2008: 
Wala wender et al.l l2008). suggesting that Bl-E is not part of an 
age gradient. 

These observations suggest that Bl-E is influenced by a 
strong, localized magnetic field. Further tests of this scenario 
would be possible by observing the entire velocity field for B 1 -E. 
Instead, we could only estimate the bulk velocity field from the 
relative motions of the individual substructures (from line cen- 
troids; see Section l4~2l ). cr^o - 0.42 km s _1 . Additionally, we do 
not have direct measurements of the gas temperature for the Bl- 
E clump. If we use the mean Bl-E dust temperature of 14 K as 
the lower limit gas temperature, c s > 0.22 km s and the veloc- 
ity field within Bl-E may be supersonic by < 2c s . These motions 
can arise in either subcritical or supercritical clouds, though the 
latter is generally associated with supersonic velocities. 

Direct measurements of the magnetic field would also be 
useful to constrain its role in Bl-E. Unfortu nately, very lit- 
tle ma gnetic field information is available here. Goodma n et al.l 
( 1990) measured optical polarization towards all of Perseus, and 
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for Bl-E, they found strong, ordered polarization vectors reach- 
ing ~ 9% polarized light, a factor of two higher than any other 
dense region in Perseus. While these polarization observations 
may suggest a strong magnetic field is associated with the Bl- 
E region, these observations are dependent on the field viewing 
angle. Also, these data were obtained at very low resolution re- 
sulting in only three vectors coinciding with Bl-E. 

A strong magnetic field does not necessarily relegate turbu- 
lent comp ression or gravity to min or roles in the evolution of B 1 - 
E. Indeed. iNakamura & Lil J2008) suggest that core forming re- 
gions exhibit different phases where magnetic fields, gravity, and 
turbulence are sometimes dominant and sometimes secondary. 
For example, they suggest that core forming regions begin with 
strong magnetic fields and strong turbulence, such that gravity is 
a secondary effect. After the turbulent energy dissipates, gravi- 
tational collapse can proceed via magnetic field lines. 

5. Conclusions 

With recent Herschel dust continuum observations from the 
Herschel Gould Belt Survey, we identified substructure in the 
Perseus Bl-E region for the first time. With our Herschel data, 
we determined the temperature and column density across the re- 
gion. We selected the nine highest column density substructures 
for complementary observations with the GBT. We summarize 
our main conclusions as follows: 

1. Bl-E contains a loose collection of roughly a dozen promi- 
nent substructures. This morphology is atypical of most star 
forming regions, which produce dense clusters and organize 
material along filaments. Such a loose collection of cores can 
be produced in magnetically subcritical sim ulations (e.g., 
see iLi & Nakamurall2002l: iBasu et al.ir2009al) influenced by 
a strong magnetic field. Furthermore, a strong magnetic field 
will delay the onset of star formation, which may explain the 
age discrepancy between Bl-E and the other nearby clumps 
in Perseus, IC 348 and NGC 1333. 

2. The Bl-E clump as a whole is gravitationally bound with 
an estimated virial parameter of a ~ 0.9. Therefore, Bl-E 
is likely to form dense cores and stars. Assuming T = 14 
K and n — 5 x 10 3 cm 3 , we find a Jeans length within the 
entire Bl-E clump of ~ 0.35 pc, accurate within a factor of 
a few. We measure a median nearest neighbour separation 
of ~ 0.15 pc for our nine substructures (~ 0.17 pc exclud- 
ing B1-E6 and B1-E7), for an inclination of 60°. Thus, the 
substructures have a median minimum separation that is less 
than the Jeans length by a factor of two. This smaller length 
scale could indicate that Bl-E contracted after the observed 
substructures formed. 

3. Several of the Bl-E substructures have masses on order 
of the Jeans mass (~ 0.9 M ), assuming T = 10 K and 
n = 1 x 10 5 cm 4 . Nevertheless, most substructures have 
large virial parameters (a » 2) indicating that they are 
gravitationally unbound. The large virial parameter suggests 
that non-thermal motions are significant and that the Jeans 
mass does not well represent the critical mass scale for these 
substructures. B1-E2 is the only substructure with a small 
virial parameter (1 < a < 3) and, due to its narrow line emis- 
sion, may be gravitationally bound. Furthermore, NH3 (2,2) 
and HC5N (9-8) were only detected towards B1-E2. Thus, 
B1-E2 is an excellent candidate for a core precursor. 

4. The Bl-E substructures have a substantial centroid veloc- 
ity dispersion (~ 0.42 km s" 1 ) resulting in an interaction 
timescale of ~ 1 Myr, assuming an average separation of 



~ 0.46 pc. This timescale is well within the lifetime of pro- 
tostellar cores. Thus, competitive accretion may play a sig- 
nificant role in the evolution of structure in Bl-E. 

The Bl-E region appears to be an excellent candidate for 
future core formation. Further studies, however, are necessary 
to explore the dust and dynamics of the bulk gas in Bl-E. 
For example, additional submillimetre continuum data along the 
Rayleigh-Jeans tail (e.g., with SCUBA-2) will improve SED- 
fitting and more accurately probe the dust emissivity index, 
/?, size, and mass. To date, our Herschel observations of Bl- 
E are the only well detected continuum data for the region. 
Additionally, our GBT observations are the only high resolution 
kinematic data for Bl-E, and we observed only the nine highest 
column density substructures seen with Herschel. These data il- 
lustrate the potential of B 1 -E to be a core forming region. Future 
observations of both the magnetic field strength and the turbulent 
velocities are necessary to probe further how Bl-E will evolve. 
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